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ABSTRACT 

CXOM31 J004252.030+413107.87 is one of the brightest X-ray sources within the D 25 region of M31, 
and associated with a globular cluster (GC) known as B135; we therefore call this X-ray source XB135. 
XB135 is a low mass X-ray binary (LMXB) that apparently exhibited hard state characteristics at 
0.3-10 keV luminosities 4-6xlO 38 erg s -1 , and the hard state is only observed below ~10% Eddington. 

If true, the accretor would be a high mass black hole (BH) (>50 M 0 ); such a BH may be formed from 
direct collapse of a metal-poor, high mass star, and the very low metalicity of B135 (0.015 Zq) makes 
such a scenario plausible. We have obtained new XMM-Newton and Chandra HRC observations to 
shed light on the nature of this object. We find from the HRC observation that XB135 is a single 
point source located close to the center of B135. The new XMM-Newton spectrum is consistent with 
a rapidly spinning ~10-20 Mq BH in the steep power law or thermal dominant state, but inconsistent 
with the hard state that we previously assumed. We cannot formally reject three component emission 
models that have been associated with high luminosity neutron star LMXBs (known as Z-sources); 
however, we prefer a BH accretor. We note that deeper observation of XB135 could discriminate 
against a neutron star accretor. 

Subject headings: x-rays: general — x-rays: binaries — globular clusters: general — globular clusters: 
individual 


1. INTRODUCTION 

One of the brightest X-ray objects in M31 is coincident 
with the globular cluster (GC) known as B135, follow¬ 
ing the Revised Bologna Catalog (RBC) v 3.4 (Galleti 
et al., 2004, 2006, 2009). We refer to this X-ray source 
as XB135. In our ^13 year Chandra survey, it exhibited 
0.3-10 keV luminosities ~4-6xl0 38 erg s _1 , assuming 
power law emissi on; the best fit photo n index (T) was 
~1.8 throughout (|Barnard et al.ll2014af) . The high lumi¬ 
nosity and hard spectrum of XB135 led us to propose 
that the accretor was a black hole (BH) with mass ~50 
Mq, which i s made possi ble by t he very low metalicity of 
XB135 ('see lBarnard et aHl2013aL and references within). 
We note that persistently bright BH binaries in GCs are 
consistent with population synt hesis theory for tidal c ap- 
ture of a main sequence star (|Kalogera et all I2004L al¬ 
though the companion might be disrupted), or with the 
formation of an ultracompa ct b inary with a de generate 
donor from a triple system (jlvanova et al.ll20Tot h 

We were awarded a 120 ks XMM-Newton observation 
of XB135, along with a 5 ks Chandra HRC observation. 
The XMM-Newton observation was designed to test the 
possibility that XB135 was formed in a low-metalicity 
environment. The purpose of the HRC observation was 
to confirm that XB135 consisted of a single point source 
that is consistent with being located in B135; XB135 had 
previously been observed at high off-axis angles, making 
its position uncertain. In this work, we fit the XMM- 


Newton spectra of XB135 with neutron star (NS) and 
BH XB emission models; additionally, we compare our 
XB135 spectrum with the XMM-Newton pn spectrum 
from the bright NS XB LMC X-2. We obtain an accurate 
position for XB135 from our Chandra HRC data. 

2. OBSERVATIONS AND DATA ANALYSIS 
2 . 1 . XMM-Newton analysis 

XB135 was observed by XMM-Newton on 2012 June 
26-27 (Obs ID 0690600401, PI R. Barnard). We analyzed 
the data using SAS version 13.0. A substantial portion 
of the observation was contaminated by background flar¬ 
ing. We screened out high background intervals in the 
pn instrument by creating a liglitcurve filtered with the 
expression “(PATTERN==0)&&(FLAG==0)&&(PI in 
[10000:12000])” and binned to 400 s; we rejected inter¬ 
vals with 10-12 keV rates >0.4 count s _1 . This resulted 
in 75 ks of good time. 

Circular source and background regions on the pn im¬ 
age were chosen, and optimized by the software. Spectra 
were extracted from these regions using the expression 
“(PATTERN<=4)&&(FLAG==0)” and the good time 
interval. We did not use any data from the MOS cameras 
because they suffered significant pile up. In our proposal 
for this observation, we made a case for using the RGS 
spectra to help constrain the metalicity. However since 
much of the observation was lost to flaring, the good time 
RGS spectra were not of sufficient quality to constrain 
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the fits. 

We estimated our uncertainties on spectral fits as fol¬ 
lows. We first obtained the best fits to the pn spectra. 
We then simulated 1000 sets of spectra from this best 
fit model using the MULTIFAKE command in XSPEC; 
variations in the simulated spectra were drawn from the 
statistical properties of the observed spectra. Each sim¬ 
ulated spectrum was modeled, obtaining the best fit pa¬ 
rameter values and the 0.3-10 keV unabsorbed flux (2-20 
keV flux in some cases, see below). Once all 1000 sets of 
spectra were analyzed, each parameter was sorted into 
ascending order, and the ler confidence limits were ob¬ 
tained from the 16 th and 84 th percentile. The parameter 
values appear to be Gaussian distributed. 

2 .2. Chandra analysis 

Our Chandra HRC-I observation was performed on 
2013 February 24 (Obs ID 14400, PI R. Barnard). We 
registered this observation to the B band Field 5 M31 im¬ 
age p r ovided b y the Local Group Galaxy Survey (LGGS, 
iMassev et alJ 120 06). using three X-ray bright GCs, an d 
following our procedure set out in lBarnard et all (2012). 
This involved using the iraf task IMCENTROID to find 
the center coordinates for each GC and its correspond¬ 
ing X-ray source, and using CCMAP to match the sky 
coordinates of the GCs in the Field 5 image to the X- 
ray positions in the Chandra image. We used PC-IRAF 
Revision 2.14.1. 

This procedure assumes that the X-ray sources are at 
the centers of their host GCs. Our assumption appears 
to be justified because the r.m.s. offset when aligning 
the X-ray and optical positions of 27 X-ray bright GCs 
in the M31 center w as jus t 0.11" in R.A. and 0.09" in 
Decl. (1 Barnard et al.ll2012ik or ~0.3 ACIS pixels. 

3. RESULTS 
3.1. Locating XB135 

The native HRC pixel size is 0.13175", and our image 
was binned by a factor 4 to ~0.5" pixels. The HRC-I 
observation revealed that XB135 is a single bright source, 
rather than a blend of multiple sources. Using the iraf 
tool IMCENTROID, we found the X-ray centroid with Im¬ 
position uncertainties of 0.07 image pixels in RA and 
Deck, i.e. 0.04". 

We found a total of four X-ray bright GCs in the field- 
of-view: B045, B116, B135, and B091D. We initially reg¬ 
istered the image using B045, B116, and B091D; this 
placed XB135 at the center of the GC B135, but the 
uncertainty in registration was negligible, perhaps due 
to our using the minimum number of sources required 
for the transformation. This resulted in a position of 
00:42:51.985 +41:31:08.32, 0.08" from the GC center in 
the Field 5 image. 

Additionally, we tried re-registering the image includ¬ 
ing XB135, to estimate the uncertainty in registration 
with an extra source; the resulting r.m.s. deviations were 
0.03" in R.A. and 0.003" in Deck The position obtained 
from this solution was 00:42:51.982 +41:31:08.31, 0.04" 
from the Field 5 GC center; the combined ler uncertain¬ 
ties (centroid position and registration) were 0.05" in RA 
and 0.04" in Deck 

For either solution, XB135 is confined to the central 
region of the GC at the 3er level; this is consistent with 



Figure 1. B band i mage of the GC B 135, from Field 5 of the 
LGGS survey of M31 (IMassev et al.|[20060 , overlaid with an ellipse 
that shows the 3cr position uncertainty for XB135, using all 4 GCs 
for registration. The uncertainties from the 3 GC solution are 
entirely contained within this ellipse. 

the expectation that the at least some BHs mi grate to 
the c enter because of mas s segregation (fMorschc r et all 
I2013T) . possibly all of them (lSpitzerll969lh FigureQ]shows 
t he GC B135 in the Field 5 B band image provided by 
IMassev et akl (120061k overlayed with an ellipse represent¬ 
ing the 3<r position uncertainties for XB135 using all 4 
X-ray bright GCs; the 3<r position uncertainties obtained 
from the 3 GC solution are entirely contained within this 
ellipse. 

We have confined a BHC to the center of a GC for 
the first time. We therefore calculated the probability 
that XB135 is an unassociated bright X-ray source that 
happens to be within 0.2" of the GC center (~10% of 
the GC radius). The XMM-Newton pipeline found 8 
sources with fluxes >1.4xl0 -12 erg s -1 in our XMM- 
Newton observation, equivalent to luminosities > 10 38 
erg s” 1 at the distance of M31. The EPIC pn field of view 
may be approximated by a circle with 14' radius, hence 
the number density of bright X-ray sources is ~4xl0 -6 
per square arc second. The probability of finding one of 
these bright X-ray sources coincidentally within 0.2" of 
the center of any_of the 80 confirmed GCs observed in the 
region (Peacock et al.ll2010f ) is ~4xl0 -5 . We therefore 
conclude that XB135 is indeed associated with the GC 
B135. We infer from this that the X-ray sources in B045, 
B116, and B091D are also located at the centers of their 
respective GCs, as suggested by our previous work. 

3.2. Lightcurves 

We present short-term and long-term lightcurves of 
XB135 in Figures [2] and [3] respectively. All lu- 
minosities assume a dist ance to M31 of 780 kpc 
ijStanek fe Garnavichlll99§ k 

Figure [2] shows the 0.3-10 keV pn intensity lightcurve 
of XB135 and the background for our 120 ks XMM- 
Newton observation. The lightcurve is binned to 400 
s and periods of high background flaring are removed. 
XB135 exhibited some variability that is not apparent 
in the background lightcurve; the r.m.s. variability is 
10.1+0.6% on time-scales longer than 400 s. The hard¬ 
ness (2.0-10 keV intensity divided by 0.3-2.0 keV in¬ 
tensity) was consistent with being constant (y 2 /dof = 
253/216), from which we infer that this variation is 
energy-independent and likely due to the stochastic vari- 
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Start Time 16104 6:33:03:182 Stop Time 16105 11:33:03:182 

Figure 2. Source (black) and background (gray) 0.3-10 keV pn 
lightcurves for XB135 from our 120 ks, 2012 June 26-27 XMM- 
Newton observation. Periods of intense background flares were 
removed. XB135 appears to have been somewhat variable, since 
it exhibited variation not seen in the background lightcurve; the 
r.m.s. variability was 10.1=b0.6%. 

at ions in accretion rate. 

Figure [3] shows a 0.3-10 keV, unabsorbed luminosity 
lightcurve of XB135, using Chandra/ACIS and XMM- 
Newton observations, represented by circles and crosses 
respectively: the ear li er XMM-Newt on observation was 
analyzed in IBarnard. Clark fc Kolb! (j2008T l. The Chan¬ 
dra ACIS lumin osities were derived from absorbed power 
law models (| Barnard et all I20121) . while the XMM- 
Newton luminosities were derived from two-component 
models. The XMM-Newton observations appear to have 
somewhat lower luminosities than the ACIS observa¬ 
tions; it is possible that the ACIS luminosities are over¬ 
estimated. However, the closest ACIS observation to our 
120 ks XMM-Newton observation occurred 15 days ear¬ 
lier, and has a luminosity that is consistent within 3cr. 
Furthermore, the variation between the archival XMM- 
Newton observation and the nearest neighbor (a fac¬ 
tor ~1.4 over 73 days) is consistent with variab ilit y we 
have s een in bright X-ray binaries (|Barnard et al.l [2012L 
l2014al) . 

3.3. Simple spectral modeling 

Since our XMM-Newton lightcurve exhibited substan¬ 
tial variability, we accumulated spectra using only the 
persistent emission (0.3-10 keV intensity <1 count s _1 ); 
this spectrum contained 53500 counts, with 2% coming 
from the background. 

We rebinned the spectrum to include at least 20 counts 
per bin. An absorbed power law fit to the spectrum 
(tbabs*po) was rejected (y 2 /dof = 1273/1033, good fit 
probability, g.f.p., = 4xl0 -7 ). An absorbed disk black- 
body fit (tbabs*diskbb) was unsuccessful also (x 2 /dof 
= 1185/1033, g.f.p. = 7xl0~ 4 ). 

Fitting a disk blackbody + power law yielded ab¬ 
sorption (IVh) = 2.7±0.3xl0 21 atom cm -2 , kTbBB = 
2.15±0.08 keV, and F = 2.2+0.3; all uncertainties are 
quoted at the la level; y 2 /dof = 984/1031 (g.f.p. = 
0.85). The 0.3-10 keV luminosity was 4.1±0.3xl0 38 erg 



Figure 3. Long term unabsorbed 0.3-10 keV luminosity 
lightcurve for XB135. Circles represent Chandra ACIS observa¬ 
tions, while crosses represent XMM-Newton observations. 

s _1 . We present this fit to the unfolded pn spectrum in 
Figure |U The 2-20 keV disk contribution to the flux 
was 69±11% (la), meaning that XB135 was not in a 
hard state (c.f. Remillard & McClintock, 2006). Instead, 
our spectrum of XB135 was consistent with either the 
thermally dominated state or the steep power law (SPL) 
state. The inner disk temperature is substantially higher 
than for most Galactic BH XBs; however GRS 1915+105 
has exhibited similar spectra, because the high BH spin 
results in a smaller I SCO and increased energy output 
(jRemillard fe McClintockll2006D . 

We note that the SPL state is only observed during 
outbursts of transient BH XBs in our own galaxy. How¬ 
ever, the two persistently bright, dynamically-confirmed 
black hole + Wolf-Rayet binaries IC10 X-l and NGC300 
X-l both exhibit emisson spectra consistent with the 
SP L state, as does the persist ent BH XB LMC X-l; in 
IBarnard. Clark fe Kolbl (120081) we suggested that such 
spectra are emitted by BH XBs with stable coronae, sim¬ 
ilar to the spectra observed from persistently bright NS 
XBs. 

One possible difference between the spectrum of XB135 
and the SPL state is that the spectrum of XB135 exhibits 
significant curvature at energies above ~2 keV. In this 
way it is similar to the high quality spectra observed 
from ultra-luminous X -ray sources (ULXs) studied by 
IGladstone et al.1 (|2009l) . Fitting our XB135 spectrum 
with a broken power law model (bknpower in XSPEC) 
had r change from 1.53+0.02 to 2.52+0.09 at 4.64+0.13 
keV; these parameters are consiste nt with several ULXs 
analyzed bv I Gladstone et ahl (I2009D . They proposed that 
ULXs exist in a special ultra-luminous state, and this 
may be true because the luminosities of the ULXs are 
considerably higher than that of XB135. However, it is 
possible that ULXs have similar spectra to other BH XBs 
(disk blackbody + power law) and the curved spectrum 
is caused by a dominant thermal comp onen t. 

We note that some authors including Gl adstone et al.1 
(120091 ) consider fits to be unphysical when the low en¬ 
ergy spectrum is dominated by the Comptonized compo¬ 
nent. However, spectral analysis of ultra-luminous X-ray 
sources in NGC253 and the confirmed BH + Wolf-Rayet 
binary IC10 X-l with XMM-Newton strongly suggests 
that this soft excess is real; fitting these spectra with the 
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Unfolded Spectrum 



Figure 4. Unfolded pn spectrum from our XMM-Newton observa¬ 
tion of XB135, with the best fit absorbed disk blackbody + power 
law model. Dotted and dashed lines represent the disk blackbody 
and blackbody components respectively. The spectrum is multi¬ 
plied by channel energy, in order to indicate the power at each 
energy. The data recieved further binning in this image for the 
sake of clarity. 


SIMP L convolut i on model for Comptonization created 
bv iSteiner et al.1 (120091) results in rejection by NGC253 
ULX3 and IC10 X-l, and results in Ay 2 = +12 over the 
two component model fits with the same degrees of free¬ 
dom for NGC253 ULX1 and NGC253 ULX2 dBarnardl 

I20l0h. 


Barnardl (|2010T ) proposed that this soft excess is due 


to an extended corona, and the evolution exhibited by 
IC10 X-l during “eclipse” provid ed very stron g evidence 
for a corona radius > 10 11 cm dBarnard et al.l [2014cT) . 
consistent with with the empirical relation between 1— 
30 keV luminosity and coro na rad iu s observed in G alac- 
tic high inclination XBs (IChurch fe Balucihska-Churchl 
12004 81. A coplanar, extended corona may feed from 
cool photons from the outer disk while heating the 
disk surface a nd stimulating production of cool photons 
(lHaardt fe Maraschilll993t ). 


3.4. BH emission models 

We decided to estimate the BH parameters for XB135 
(assuming a BH accretor) by adopting a model that has 
been used to describe Galactic BH bin aries. Stu d ies o f 
the metalicity in the M31 bulg e dJablonka &; Saraiedinil 
2005) and part of the halo (|Durrell. Harris fe Pritchetl 
20011 ) show that the metalicity varies over 2 orders of 
magnitude in each region. While it is unlikely that the 
host GC B135 (Z ~ 0.015 Zq) harbors much material, we 
obtained best fit spectra for models where the absorber 
metalicity Z = 0.1-1.0 Zq as well as Z = 0.015. To do 
this, we altered the metalicity with the ABUND command 
in XSPEC. 

The XSPEC model used was 
WABS*TBABS*SIMPL*KERRBB, described in 
the following paragraphs. For absorption, we included 
7xl0 20 atom cnr 2 at Solar metalicity to account 
for Galactic absorption (WABS); we also included an 
absorber with variable metalicity (TBABS, with the 
abundance controlled by the ABUND command). 

We modelled the emission as SIMPL*KERRBB, de¬ 
scribing the partially Comptonized emission from a disk 


around a spinning BH (following ISteiner et al.l 12009D . 
The variables for the emission model are: the fraction 
of disk emission that is Comptonized (fc), the BH spin 
parameter (a), the disk inclination (i), the BH mass 
(Mbh), and the accretion rate (M, normalized to 10 18 
g s -1 ). The KERRBB model also requires a distance 
(assumed to be 780 kpc); the spectral hardening factor, 
self-irradiation flag and limb darkening flag of the KER¬ 
RBB model were kept at the default settings. For the 
Comptonized emission, we assumed T = 2.4, as is typi¬ 
cal for the steep power law state. In the end, the results 
proved rather insensitive to T. 

Decreasing the M31 metalicity resulted in hotter spec¬ 
tra, translating as lower BH masses. The solar metalicity 
fit yielded the highest BH Mass: 21.8±0.8 Mq , with to¬ 
tal Ah = 1.60±0.04xl0 21 atom cm” 2 , f c = 0.57+0.03, 
a = 0.998+0.005, i = 65.7±0.7°, M = 8.93±0.05xl0 17 
g s -1 , and \ 2 /dof = 1008/1030. Assuming an absorbed 
metalicity of 0.015 outside our Galaxy yielded a BH 
mass of 16±3 Mq, with A // 731 = 4.9±0.4xl0 21 atom 
cm 2 , f c = 0.34+0.10, a = 0.95+0.02, i = 78±3°, M = 
2.5±0.5xl0 18 g s _1 and y 2 /dof = 977/1030. All uncer¬ 
tainties are quoted at the ler level. 

We caution that these results should not be treated 
as accurate, since we make assumptions in our modeling 
that may not apply. As such, we find that the mod¬ 
els suggest a ~10-20 Mq BH with spin ~1, somewhat 
like the Galactic BH LMXB GRS 1915+105 which has 
a spin >0.98 dMcClintock et al.l l2006f ). They only con¬ 
sidered spectra with luminosities <30% Eddington when 
estimating the spin, to ensure that the disk is thin. Our 
Comptonized Kerr BH models yielded 0.3-10 keV lumi¬ 
nosities ~3xl0 38 erg s _1 , ~12-24% Eddington for BH 
masses ~ 10-20 Mq, whi ch is in the range considered by 
IMcClintock et al.1 (|2006f) : our disk blackbody + power 
law fit yielded a luminosity ~4xl0 38 , which is 16-32% 
Eddington for a 10-20 Mq BH. Unlike GRS 1915+105, 
XB135 appears to have a strong Comptonized compo¬ 
nent, which may co mplicate our interpretati on. 

We also note that IMcClintock et al.l (120061) identified a 
weakness in the KERRBB model: the spectral hardening 
parameter (color temperature divided by effective tem¬ 
perature) is assumed to be a constant. Instead, they used 
a tabular model KERRBB2 where the hardening is com¬ 
puted. However, the hardening factor for the range of 
parameters that we are interested in is likely to be ~1.7, 
the value assumed in our model (J. Steiner, private com¬ 
munication). Hence the differences between KERRBB and 
KERRBB2 models are not expected to be great for XB135. 

3.5. Comparison with Galactic NS binaries 

Gal actic NS XBs exhibit a wide range of behaviors; re¬ 
cently [LnTeFal] (l2007Ll2009ll2012t) parameterised the full 
gamut by examining thousands of RXTE spectra from 
persistent and transient XBs exhibiting all spectral states 
using the same disk blackbody + blackbody emission 
model. They found that their model described all NS 
XB spectra well, with two exceptions: low accretion rate 
state s where the emission dominated by Comptoniza¬ 
tion (|Lin et al.l 120071I2009D : and the so-called Horizon¬ 
tal Bran ches of the highest luminosity NS XBs called Z- 
sources (|Hasinger fc van der Klisl 119891 ) where a Comp¬ 
tonized component is required in addition to the disk 
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black b ody and blackbo dy components dLin et al.l I2009L 
I2012D . ILin et all (120101) successfully applied their double 
thermal model to Beppo-SAX and Suzaku spectra of a 
persistent NS XB, with energy ranges 1.0-150 keV and 
1.2-40 keV respectively; hence the apparent success of 
the double thermal model is not confined to the RXTE 
band. We note that this approach is not necessarily real¬ 
istic, but the approach is valuable because it allows us to 
assess the full gamut of NS LMXB behavior in a single 
parameter space. 

We have identified 50 BHCs in M31 by comparing disk 
blackbody + blackbody fits to their spectra with the 
range of parameters exhibited by Galactic NS XBs (see 
iBarnard et al.ll2014bl and references within). For each 
BHC, we calculated the probability that the disk black¬ 
body temperature, blackbody temperature, and the disk 
blackbody contribution were consistent with the NS XB 
parameter space (Pdbb, +bb and Pf respectively). Mul¬ 
tiplying these probabilities together yielded the proba¬ 
bility that the BHC was consistent with being a NS XB 
(.Pns)- Since the only Galactic NS XBs to exhibit lumi¬ 
nosities comparable with XB135 are Z-sources, we com¬ 
pared our spectrum with the Z-source states studied by 
ILin et all (1200911201211 : they found that A;Tdbb >1-3 keV, 
kTBB >2.2 keV, and /d >0.66. 

3.5.1. Comparison with disk blackbody + blackbody emission 

models 

Lin et al. (2009; 2012) found that the major¬ 
ity of spectra were well described by disk blackbody 
+ blackbody emission models. We created mod¬ 
els with Z = 0.015-1. OZq using the XSPEC model 
WABS*TBABS*(CFLUX*DISKBB+CFLUX*BB); the 
fluxes were obtained in the 2-20 keV range for better 
comparison with the results of Lin et al. (2009; 2012). 

As the metalicity of the M31 absorber increased, 
fcTbBB, kTBB, and /b systematically decreased; fcXbBB 
fell from 1.15+0.10 keV to 0.87+0.05 keV, kT BB from 
1.84±0.15 keV to 1.57±0.05 keV, and / D from 0.41±0.08 
to 0.26±0.03. Pns ranged from 10 -5 ' 1 to 10” 93 . 

3.5.2. Comparison with disk blackbody + blackbody + power 

law emission models 

Lin et al. (2009; 2012) found that the Horizon¬ 
tal branch states of XTE J1701462 and GX 17+2 re¬ 
quired a Comptonized component, and found that the 
combined flux contribution of the disk blackbody and 
Comptonized component (/dbb+pl) was >76% of the 
total. Hence we also made a series of models with 
WABS*TBABS*(CFLUX*DISKBB+CFLUX*BB+ 
CFLUX*PO); we fixed the photon index of the Comp¬ 
tonized component to 2.4, as before. 

For these models /cTdbb and kTBB systematically 
increased with increasing metalicity: ATdbb from 
1.14+0.10 keV to 1.4+0.4 keV, and kT B b from 1.83+0.15 
keV to 2.4+0.8 keV, with /dbb+pl rising from 0.35+0.05 
to 0.6+0.3. We find that for Z >0.4, Pns >0.01, and 
XB135 is consistent with NS XB spectra fitted with this 
emission model at the 3er level. However, we note that 
the uncertainties in all parameters are significantly larger 
for this three-component emission model than for the 
previous two-component model; hence deeper observa¬ 
tions of XB135, or with a broader energy range, could 
potentially separate XB135 from Galactic Z-sources for 


the three component model as well as the two component 
model. 

3.6. Comparison with an XMM-Newton observation of 
the Z-source LMC X-2 

LMC X-2 is a NS XB residing in the Large Magellanic 
Cloud which exhibits spectral and timing behavior rem¬ 
ini scen t of a Sco- like Z-source (ISmale fc Kuulkeri|200Ct 
ISmale. Homan, fe Kuulkersl 120031) . It is also extremely 

bright for a NS XB, ~0.5-2 Leap where +Edd is the Ed¬ 
dington li mit. It has a known co unterpart, and a mass 
ratio <0.4 (|Cornelisse et al.lfeOOTl) . 

Unlike the Galactic Z-sources, LMC X-2 is suf- 
ficiently faint for o bservation with XMM-Newton; 
ILavagetto et al.l (12008D present their analysis of the 2003 
April 21 XMM-Newton observation of LMC X-2, which 
was conducted in Small Window mode with the Medium 
filter, allowing pn spectra that were free from pile-up. 
They used SAS version 7.0 and XSPEC version 11.3.2ad 
to analyze the pn and RGS spectra. After establish¬ 
ing that the background was stable, and tha t the s pec- 
trum varied little over time, ILavagetto et al.l (|2008l) ex¬ 
tracted source and background spectra from the whole 
observation. They found a large discrepancy between 
the pn and RGS spectra for energies below 1.3 keV, 
and attribut ed this to poor pn cali bration for bright X- 
ray sources dBoirin fe Parmanl2003D : they therefore only 
u sed pn data abov e 1.3 k eV. 

ILavagetto et al.l (120081) fitted the pn and RGS spec¬ 
tra with various models, including the disk blackbody 
+ blackbody model adopted by Lin et al. (2007, 2009, 
2012). Their best fits yielded a disk blackbody tempera¬ 
ture of 0.815+0.002 keV, and a blackbody temperature of 
1.543+0.009 keV, which are significantly cooler than the 
parameters obtained from the Z-source spectra of XTE 
J1701—462 and GX 17+2 (ILin et al]l2009t l20ll) and re¬ 
markably similar to the parameters that we obtained for 
XB135. We therefore decided to investigate the XMM- 
Newton observation of LMC X-2 ourselves. 

Using SAS version 13.0 and XSPEC version 12.8.Ip, 
we modeled the pn spectrum above 1.3 keV using a 
disk blackbody + blackbody model (XSPEC model 
WABS*(DISKBB+BB) for consistency with Lin et al. 
(2007, 2009, 2012)). The spectrum contained ~750,000 
counts, with the source contributing ~96%. We binned 
the spectra to a minimum of 20 counts per bin. The ab¬ 
sorption was unconstrained, so we adopted the bes t fit 
column density obtained by ILavagetto et al.l (I2008D for 
the disk blackbody + blackbody model, 4xl0 20 atom 

_ 9 

cm . 

Our best fit disk blackbody and blackbody tempera¬ 
tures for LMC X-2 were 1.36+0.04 keV and 2.42+0.08 
keV respectively with y 2 /do f = 1812/1737, very d iffer- 
ent to the results obtained bv lLavagetto et all (I2008D and 
consiste nt with the paramete rs observed in Galactic Z- 
sources (ILin et al.ll2009l I2012D . When t he temperatures 
were frozen at the values obtained by ILavagetto et all 
M) , y 2 /dof = 3002/1739. 

We present our best fit disk blackbody + blackbody 
fits to 0.3-10 keV and 1.3-10 keV pn spectra for LMC 
X-2 in the left and right panels of Figure [5] respectively. 
We also present residuals. The 0.3-10 keV fit yielded 
unaccepatable results, with strong systematic variations 
in the residuals. The null hypothesis probability for the 
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Figure 5. Best fit disk blackbody + blackbody models for the 0.3—10 keV pn spectrum (left) and 1.3-10 keV pn spectrum (right) from 
the Z-source LMC X-2. The dotted and dashed lines represent the disk blackbody and blackbody components respectively. The lower 
panel shows the x 2 residuals. For the left panel, we see considerable oscillation in the residuals, indicating a bad fit from the 0.3-10 keV 
spectrum. For the 1.3-10 keV fit, there are no systematic residuals, but some channels show large deviations that are probably due to 
uncertainties in calibration in this 750,000 count spectrum. 


1.3-10 keV fit (i.e. that differences between model and 
data are due to random variations in the spectrum) is 
0.10, making this an acceptable fit even though \ 2 ~5- 
10 for some channels; these are likely due to uncertainties 
in calibration, which are only significant for very bright 
sources. 

The fact that our new parameters for LMC X-2 are 
consistent with the Galactic Z-sourc e s sugg ests that the 
results obtained by lLavagetto et ahl (120081) were proba¬ 
bly affected by calibration errors, especially since they 
reported the discrepency between RGS and pn fits be¬ 
low 1.3 keV vet chos e to fit the RGS spectra. We note 
that Lavagetto et al.l ( 200 8) reported >5400 degrees of 
freedom in their spectra, and most of these come from 
low-weighted RGS bins; hence it is possible that their 
RGS+pn fit would be rejected if they had only consid¬ 
ered the pn data. 

Fitting an absorbed power law model to the 1.3-10 
keV pn spectrum of LMC X-2 allows us to compare the 
general shape of its spectrum with that of XB135. We 
were unable to constra in the absorpti o n due to the 1.3 
keV cut-off; however, lLavagetto et al.1 (120081) found Nu 
to be ^4-lOxlO 20 atom cm -2 , depending on the emis¬ 
sion model. For Ah = l.OxlO 21 , T = 1.373±0.003, with 
y 2 /dof = 4546/1739, and lower values of Ah yielded 
harder T; such a spect rum is significantly harder than 
any BH spectra (|R.emillard fc McClintockI 12006 ). Sim¬ 
ilarly, fitting XB135 with the best fit temperatures of 
LMC X-2 results in % 2 /dof = 999/1033, whereas the 
equivalent free fit yielded y 2 /dof = 973/1031; F-testing 
shows that the probability that this improvement is due 
to random variations is 1.2xl0~ 6 , and the spectra of 
XB135 and LMC X-2 differ at a >4.8er level. The spec¬ 
trum of XB135 is significantly softer than that of LMC 
X-2 when observed with the same instrument, strength¬ 
ening our case for XB135 containing a BH accretor. 

3.7. Comparison with MSI XBs 

While it is instructive to compare the X-ray popula¬ 
tions of M31 and the Milky Way, this does not provide 
the full picture because the X-ray populations of these 
two galaxies are very different. 


iStiele et~aL] (120111) found 1897 X-ray sources within 
the M31 D 25 region at luminosities > 10 35 erg s _1 ; 
around 60% of these sources were either class ified only 
as Hard, or unclassified, and IStiele et al.l (| 2 Qlll ) expected 
these to be AGN. However, the Chandra survey of the 
M31 center (within ~20' or 4 kpc of M31*) revealed 
over 200 X-ray sources with luminosities >10 35 erg s " 1 
that were missed by XMM-Newton, th anks to Ch an¬ 
dra’s l ower background and superior PSF (Barnard et al. 
I2014al) . Furthermore, we identified ~200 new X-ray bi¬ 
naries, and found strong evidence that the population of 
X-ray sources consistent with AGN contains a number 
of unidentified XBs. Hence there are more X-ray sources 
an d fewer AGN within the M31 D 25 region than expected 
bv IStiele et al.1 (1 201 lh. 

In iBarnard et al.l (|20141j ). we present analysis of 50 
BHCs in M31, adding 13 transients and 2 GC XBs stud¬ 
ied with XMM-Newton to our existing sample. By com¬ 
paring the peristent and transient XB population within 
6 ' of M31*, we found that >40% of XBs >10 35 erg 
s _1 are expected to contain BH primaries. Remarkably, 
BHC XBs represent >90% of those XBs with luminosi¬ 
ties >10 38 erg s _1 ; this is very different to the Milky 
Way, where the majority of XBs >10 38 erg s _1 contain 
NS primaries. 

Of the 15 BHCs associated with GCs, 12 exhibit disk 
blackbody + blackbody spectra th at differ from NS spec¬ 
tra by >5cr ()Barnard et al.ll2014bl ). While ~30 M31 GCs 
contain X-ray sources that have exceeded 10 37 erg s _1 at 
some point in the last 13 years, only one Galactic GC 
has ever been observed to exceed this luminosity. Hence 
the lack of confirmed BH XBs in Galactic GCs does not 
hinder our classification of XB 135 as a BHC. 

We know of two M31 X-ray sources that exceed 10 38 
erg s _1 in the 0.3-10 keV band, and most likely con¬ 
tain NS accretors. First is XB158, a high inclination 
X-ray binary that exh ibits repeated intensity v ariations 
on a ~2.8 hr period (iTrudolvubov et ahl 12002 ) as well 
as a superorbital period of ~5.8 days that was recently 
revealed by daily Swift observations over 30 days (R. 
Barnard et al. 2015, ApJ submitted). Its 2002 XMM- 
Newton spectrum is consistent with disk blackbody + 














































XMM-Newton and Chandra study XB135 


7 



o LJ-.-1- 1 -1- 

0 2000 4000 

MJD - 51450 / day 



LI_i_i_i_i_ 

0 2000 4000 

MJD-51450/day 


Figure 6. Chandra ACIS 0.3-10 keV luminosity lightcurves (top) and spectral histories (bottom), for the M31 Z-source RX J0042.6+4115 
(left) and XB135 (right), presented with equal scaling. We found the best fit absorbed power law model for each observation, and examined 
the spectral variability by fitting a line of constant photon index; XB135 is consistent with no spectral variation at the 3 <j level, while the 
Z-source RX J0042.6+4115 exhibited substantial spectral variability. 


blackbody + power law fits to Galactic NS XBs; how¬ 
ever, a simple power law fit yields T 0.57±0.09, con¬ 
siderably harder than any BH spectrum (Barnard et al. 
2015). Second is RX J0042.6+4115, which exhibited tri- 
mo dal spectral evolution analagous to Galactic Z-sources 
(|Barnard et al.ll2003l) . Power law fits to our ~100 Chan¬ 
dra ACIS observations of RX J0042.6+4115 yielded a 
mean T of 0.856+0.008 (|Barnard et al.ll2014alf . although 
the spectrum varied considerably, as discussed below. 
While these power law emission models are clearly not 
realistic, they demonstrate that these two bright M31 NS 
XBs are significantly harder than XB135. 

In Figurc[H]we compare the 0.3-10 keV lightcurves and 
spectral evolution from XB135 and the M31 Z-source 
candidate RX J00421.6+41152 over 13 years’ worth of 
Chandra (see iBarnard et all l2014al for main survey 
results). Z-sources evolve d over time-scales of days to 
weeks (|Muno et al.l I2002H . hence if XB135 were a Z- 
source, we should expect substantial spectral variabil¬ 
ity. We present the 0.3-10 keV liglitcurve in the up¬ 
per panel, and fit photon index for the best fit absorbed 
power law emission model in the lower panel; the scales 
are identical. We fitted each object with the best fit 
line of constant T. For XB135, x 2 /dof = 27/12, with 
probability 7xl0~ 3 , meaning that T is consistent with 
being constant at the 3<r level. For the Z-source RX 
J0042.6+4115, the best fit line of constant T yielded 
X 2 /dof = 381/112, with a probability of 2xl0" 31 . We 
note that the ACIS observations of RX J0042.6+4115 
did result in pile-up; however, one would expect T to de¬ 
crease for a piled-up source as the intensity increased, 
yet r for RX J0042.6+4115 substantially increased, so 
the true spectral variation of RX J0042.6+115 was likely 
to be greater than observed. We present this as evi¬ 
dence in support of a BH accretor for XB135. We note 
that RX J0042.6+4115 exhibited flares during the 2002 
XMM-Newton observation similar to those observed in 
Figure 2; however, the flares observed in RX J0042.6+15 
were clearly ener gy dependent w ith higher amplitudes at 
higher energies (IBarnard et al.1 12003T) . in contrast with 
the energy-independent flares observed from XB135. 


4. DISCUSSION AND CONCLUSIONS 

We obtained new XMM-Newton and Chandra obser¬ 
vations of the BH candidate associated with the M31 
globular cluster B135. Previous results were suggestive 
of an extremely high BH mass of ~50 M 0 . The 120 ks 
XMM-Newton yielded ~50,000 net source counts from 
intervals free from background flares and source varia¬ 
tion. A high luminosity, hard state fit to such a high 
quality spectrum would have been extremely compelling 
evidence for such a massive accretor; unfortunately, such 
a hard state was ruled out. Fitting a black hole emission 
model to our XB135 spectrum suggests a black hole mass 
~10-20 M 0 , considerably smaller than the ~50 M@ BH 
previously predicted. 

Our Chandra HRC observation provided the first ob¬ 
servation of XB135 where the PSF was not greatly en¬ 
larged due to being observed at a high off-axis angle. We 
find that XB135 is a single point source, located at or 
very near the center of the GC B135 with a 3cr position 
uncertainty of 0.15" in R.A. and 0.12" in Deck 

We prefer a BH primary for reasons stated earlier in 
the paper. We expect >90% of M31 X-ray binaries with 
0.3-10 ke V luminosities >10 38 erg s _1 to contain BH 
accretors (iBarnard et al.lf2014bl) , and ~90% of those very 
bright X-ray sources are located in GCs. Also, the 0.3- 
10 keV spectrum of XB135 is significantly softer than 
the similarly luminous Z-source LMC X-2 when studied 
with the same instrument (XMM-Newton pn), as well 
as RX J0042.6+4115 and XB158. Furthermore, XB135 
exhibits less spectral variability than the M31 Z-source 
RX J0042.6+4115. Fitting the XB135 spectrum with a 
disk blackbody + power law model yields results that 
are consistent with the steep power law state when the 
uncertainties are considered; while the steep power law 
state is only exhibited by transient BH XBs in our own 
galaxy, persistent BH XBs IC10 X-l, NGC300 X-l, and 
LMC X-l all exhibit spectra consistent with the steep 
power law state. 

However, some aspects of the behavior observed from 
XB135 may contradict our BH assessment. We note 
that the best fit disk blackbody + power law param¬ 
eters for XB135 are not obviously attributable to any 
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canonical black hole state; however, the best fit disk 
blackbody + blackbody + power law model for XB135 
is consist ent with the Galactic Z-source studied by 
iLin et al ] Em [ 20 T 1 . We also note that the type of 
flaring behavior observed from XB135 (see Fig. 2) is 
common for Z-sources, but not for BH XBs; however, 
the energy-independent flaring behavior from XB135 
may be distinctly different from the energy-dependent 
flares observed from XMM-Newton observations of RX 
J0042.6+4115 (IBaxnard et a,l.ll200l . 

With the nature of the accretor in XB135 unclear, fur¬ 
ther observations are required. The uncertainties on the 
three component model would be constrained by either 
deeper observations, or by spectra with wider spectral 
ranges. Also, a better understanding of the metalicity of 
the absorbing material in M31 would also constrain the 
spectral parameters, and may aid the classification. 
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